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Abstract 

We performed Monte Carlo simulations of the two-dimensional q- 
state Potts model with q = 10, 15, and 20 to study the energy and 
magnetization cumulants in the ordered and disordered phase at the 
first-order transition point (3f By using very large systems of size 
300 X 300, 120 X 120, and 80 x 80 for q = 10, 15, and 20, respectively, 
our numerical estimates provide practically (up to unavoidable, but 
very small statistical errors) exact results which can serve as a useful 
test of recent resummed large-g' expansions for the energy cumulants by 
Bhattacharya et al. [J. Phys. I (France) 7 (1997) 81]. Up to the third 
order cumulant and down to (7 = 10 we obtain very good agreement, 
and also the higher-order estimates are found to be compatible. 
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1 Introduction 


A detailed understanding of first-order phase transitions plays an important 
role in many fields of physics [1]. In particular the hnite-size scaling behaviour 
near the transition has been a snbject of increasing interest in recent years [2], 
It is well known, and in limiting cases exactly proven [3], that thermodynamic 
observables in equilibrinm can be expanded as asymptotic power series in 
l/I/, where V is the volume of the system. The range of applicability of such 
expansions, however, is widely nnknown and was the sonrce for quite a few 
apparent inconsistencies in the recent literatnre. As has been discussed in a 
recent series of papers by Bhattacharya et al. [4-8], a precise knowledge of 
the energy cnmnlants of the coexisting phases at the transition point (3t in the 
inhnite-volume limit can help to understand and to resolve these problems. 

The explicit calculations by Bhattacharya et al. have been performed 
for the two-dimensional (2D) g-state Potts model [9, 10], which for g > 5 
is the paradigm for a system with a temperature-driven first-order phase 
transition. The advantage of this model is that many properties (transition 
point Pt [11]; internal energy [11] and magnetization [12] of the pnre phases at 
Pt, correlation length [13] and interface tension [14] at Pt) are exactly known 
and that the parameter g (the nnmber of states per spin) allows one to tnne 
the strength of the transition. Based on the Fortnin-Kasteleyn representation 
[10, 15] of the Potts model, Bhattacharya et al. analyzed the ensuing clusters 
to obtain large-g series expansions for the energy cumulants in the ordered 
phase at Pt and then applied Pade-type resummations to arrive at nnmerical 
estimates. 

Since, as with every resnmmed series expansion, it is intrinsically diffi¬ 
cult to provide reliable estimates of the size of systematic errors, we found it 
worthwhile to determine the cumulants by a completely independent method, 
namely Monte Carlo simulations. For large enough system sizes the system¬ 
atic errors are negligible (of the order exp{—L/^), where L is the linear size 
of the system and ^ the hnite correlation length), and by increasing the simn- 
lation time also the statistical errors can be made as small as desired. In the 
following we report high-statistics measnrements of the first ten energy cn- 
mulants in the ordered and disordered phase at Pt for the models with g = 10, 
15, and 20. With increasing order we observe the expected steep growth of 
the energy cumulants. The accnracy of the estimates, however, decreases 
with increasing order because the tail ends of the energy distribution become 
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more and more important. Up to the third-order cnmnlant the precision of 
the nnmerical estimates is still high (abont 1% or better), and we obtain very 
good agreement with the large-g expansions down to g = 10. At fonrth order 
the agreement between the two methods is still good within the statistical 
error bonnds of abont 5% — 15%, and also the results for the hfth and sixth 
order, the highest energy cumulants considered in Ref. [8], are compatible 
with each other. The further estimates up to the tenth order should only be 
taken as a rough indication of the order of magnitude in high orders of the 
cnmnlant expansion. In addition we also present estimates of the hrst three 
magnetization cumulants in the two phases at I3t, and compare the suscepti¬ 
bility in the ordered phase with recent low-temperature series expansions for 
q = 10. 

The remainder of the paper is organized as follows. In Sec. 2 we hrst 
briehy recall the model and some exact results. We then describe the set-up 
of our simulations and discuss the estimators used for measuring the various 
cumulants. The numerical data are presented in Sec. 3, and in Sec. 4 we 
conclude with a brief summary of the main results and a few hnal remarks. 

2 Model and simulation 

In our simulations we employed the standard dehnition of the Potts model 
partition function [9, 10] 

^ ^ = - H Si = l,...,q, (1) 

{«d (o> 

where (3 = J/ksT is the inverse temperature in natural units, F is the 
free energy, i denote the sites of a two-dimensional square lattice, {ij) are 
nearest-neighbor pairs, and 5sis is the Kronecker delta symbol. 

All simulations were performed in the canonical ensemble at (3 = (3t = 
ln(l -|- y/g), employing large square lattices of size 300 x 300, 120 x 120, and 
80 X 80 for g = 10, 15, and 20, respectively, and periodic boundary conditions. 
To stabilize the pure ordered or disordered phase we took advantage of the 
extremely small tunneling probability for large system sizes [16-19]. Since 
the tunneling proceeds via mixed phase conhgurations with two interfaces, 
the probability is proportional to exp{—2aodL), where aod is the interface 
tension between the ordered and disordered phase. For 2D Potts models 
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aod has been analytically predicted [14], 2aod = l/'Cd(A)) where ^d{Pt) is 
the exactly known [13] correlation length in the disordered phase at the 
transition point (= 10.559519..4.180954..., and 2.695502... for g = 10, 
15, and 20). By rewriting the tunneling probability as exp{—L/^d), h is 
easy to estimate that for our choice of lattice sizes, L ^ 28^d, the order of 
magnitude is about exp(—28) ~ 10“^^. It is therefore extremely probable 
that, starting from a completely ordered or disordered configuration, the 
system will stay in the ordered or disordered phase during a very long (but 
finite) simulation time, thereby allowing statistically meaningful pure-phase 
measurements of energy and magnetization cumulants. To be sure we have of 
course monitored the time evolution of our simulations and explicitly verified 
that no tunnelings occurred. In Fig. 1 we show the probability distributions 
P(e) of the energy density e = E/V in the ordered as well as in the disordered 
phase, demonstrating that the two peaks are indeed very well separated. 

The finite-size corrections in the pure phases are also expected to be of 
the order exp(—L/,^p), where the subscript p stands for the ordered (o) and 
disordered (d) phase, respectively. Since we have recently obtained strong 
numerical evidence that ^o{Pt) = ^d{Pt) [18, 19], this yields for the cho¬ 
sen lattice sizes in both phases again an order of magnitude estimate of 
exp(—28) 10“^^. 

To update the spins we employed in the ordered phase the heat-bath 
algorithm, while in the disordered phase it is more efficient to use the single¬ 
cluster algorithm [16, 17]. The observed integrated autocorrelation times 
of the energy, rint,e; and the statistics parameters are compiled in Table 1. 
Notice that for the heat-bath algorithm rint,e scales roughly with as one 
would have expected on general grounds. For the single-cluster algorithm we 
followed the usual convention and defined I^/(|C|)rf,sc single-cluster steps as 
one Monte Carlo update sweep (MCS), where (|C|)rf,sc is the average cluster 
size (cf. Table 5 below). 

The simulations with the heat-bath algorithm were performed on a CRAY 
vector computer. To estimate statistical errors we divided the runs into 
several bins and employed the standard jack-knife procedure. The single¬ 
cluster code was implemented on a T3D parallel computer by simulating 64 
time histories in parallel. This enabled us to gather an equivalent of about 
five workstation CPU years within a relatively short time. Here the error bars 
are estimated from the fluctuations between the 64 copies by using again the 
jack-knife procedure. 
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The primary observables we report here are the energy cumulants 
at the transition point f3t = ln(l + ^/q), which are dehned through the Taylor 
expansion of the scaled free energy density (3f = f3F/V around f3 = Pt-, 

- lifpiP) = -A/(A) + I:(-i)” 4”’(A)(/3 - A)7n! ■ (2) 

n=l 

While the free energy is continuous at a hrst-order phase transition, f{Pt) = 


foiPt) = fdiPt), the cumulants are discontinuous and hence different in the 
ordered and disordered phase. The hrst three cumulants coincide with the 
(central) moments, 

= e,= {E)JV, (3) 

(4) 

4"' = ^{(E-{E)p)\/V, (5) 

where Cp is the usual specihc heat in the pure phases. The higher-order cu¬ 
mulants can be expressed as non-linear combinations of the central moments 
4 "’ = ifE - {E),r),/V, e,g,, 

4'*' = (6) 

4^^' = 4"-ioWW. (7 

4®' = 4®’- i5r4®’4'‘’- ior4’>" + (s) 


and so on as listed up to the tenth order in the Appendix. While at Bt both 
Co and Cd are known exactly [11], the energy cumulants and with 
n>2 can only be related to each other via duality [10]. In particular for Co 
and Cd as well as and the duality relations read 

Co = Cd - Pticd- eo)/y/q, (9) 

= -/^?^+2(l-g)/g^'^^-3(ed-eo)/g + 6cd/(A^yg). (10) 

For the magnetization cumulants we have used slightly different deh- 
nitions in the ordered and disordered phase. In the disordered phase the 
magnetization vanishes and the magnetic susceptibility can be dehned as 
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which is certainly independent of the choice of the reference orientation 
s = 1,... ,q. In simulations employing cluster-up date algorithms the same 
quantity can also be estimated from 


Xd = q_ I (l^l)t^,sc 


1 (|C*P)d,SW 

Q (|C'|)d,SW ’ 


( 12 ) 


where \C\ denotes the size or weight of a cluster, i.e., the number of spins 
belonging to a cluster, and the subscripts “SC” and “SW” refer to averages 
over the clusters encountered in the single-cluster and Swendsen-Wang for¬ 
mulation, respectively. Alternatively one could also measure the (projected) 
spin-spin correlation function ( 7 (f), using the spin or one of the cluster repre¬ 
sentations, and derive Xd as the “integral” of g{i), Xd = ~ 1 )^) HiQii) 

[16, 17], 

In the ordered phase we have measured the maximum definition of the 
magnetization 

M^ax= iV^ax = max{iVi,iV 2 ,...,iVj, (13) 

where Ns counts the number of spins of orientation s = 1 ,..., g in a given 
configuration. A cluster estimator for the magnetization is 


Afclus — |Cjnax|) (14) 

where |Cmax| denotes the size of the largest (spanning) stochastic cluster 
in each spin configuration. The expectation values m = (Mmax)o/f^ and 
m' = (Mcius)o/C coincide, m = m', and can be directly compared with the 
exact result for g > 5, m = [(1 — a;"')/(l — x^”)], with 0 < a; < 1 defined 

by g = a; -|- 2 -|- x~^ [12]. The magnetic susceptibility in the ordered phase is 
computed as 

Xo = ({aO„ - /c (15) 

and the third-order magnetization cumulant (or central moment) is given by 

= ((Mmax - {Mraa^)of)o/V. (16) 


Notice that while the expectation values of M^ax and Mdus coincide, this is 
not the case for the susceptibilities Xo and y),, where the latter is defined by 
(15) with Mmax replaced by Mdus- The same remark applies of course also 
to the third-order cumulants and 
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3 Results 


The Monte Carlo results for Cp, Cp, and /i® at jSt are collected in Table 2 
{p = d, disordered phase) and Table 3 (p = o, ordered phase). For comparison 
we have also listed previous measurements [17, 19] on smaller lattices which 
clearly demonstrate that the above described “dynamically stabilized” pure 
phase simulations are feasible and that, as expected, the residual hnite-size 
corrections are completely covered by the statistical errors. 

The Table 2 shows the Monte Carlo results for the first three energy 
cumulants at f3t obtained from the simulations in the disordered phase. 
The hrst moments, e^, are in excellent agreement with the exact results, 
with relative statistical errors of the order 2 x 10“^. This indicates quan¬ 
titatively that systematic errors are completely under control (including 
possible problems with pseudo-random numbers). The lines denoted by 
“(large g)” show for the readers convenience the duality transformed esti¬ 
mates from the large-g expansions for the cumulants in the ordered phase 
(cp. Table 3). To be specihc, for /i® we have rewritten (10) with the help 
of (9) as = -/i® -h 2(1 - g)/g^/^ -h 3(6^ - eo)/q + and in¬ 

serted the large-g estimates of Table 3, which are taken from the most recent 
publication [8]. 

In Table 3 we show the first three energy cumulants at f3t in the ordered 
phase obtained first directly from the simulations in the ordered phase (de¬ 
noted by, e.g., (o, L x L)), and second via the duality relations (9) and (10) 
using the just described results in the disordered phase (denoted by, e.g., (d, 
L X L)) collected in Table 2. Also in the ordered phase the exactly known 
hrst moments are conhrmed with high precision, and the duality relations are 
very well satished by the two independent sets of Monte Carlo simulations, 
which further underlines the reliability of the data. As already mentioned 
above, the large-g expansion estimates are taken from Ref. [8]. (Thanks 
to longer series expansions and constantly improved analysis techniques the 
numbers given in [4-8] scatter a little bit, rehecting the current state of the 
art.) We see that in all cases the agreement between the Monte Carlo and 
large-g estimates is extremely good. The only exception is perhaps ^ for 
g = 10, but here quite naturally the systematic error of the large-g expansion 
is already relatively large. Recent analyses of low- and high-temperature se¬ 
ries expansions for the g = 10 Potts model specihc heat at on the other 
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hand, yielded much larger values of Co = 31.8(2.8) and q = 33(3) [20]. 

For the 2L x L and 2L x 2L lattices we have also computed higher- 
order energy moments and the resulting cumulants. Our results up to the 
eighth order are collected in Table 4. With increasing order these observables 
become very sensitive to the tail ends of the energy distribution and the 
statistical accuracy deteriorates quite rapidly. Due to cancellation effects 
this decrease in accuracy is much more pronounced for the cumulants than 
for the (central) moments. For the statistical errors are about 5% —15%, 
and here we still hnd quantitative agreement of with the resummed large- 
q expansions, whose estimated systematic errors are of the same order. Also 
for and the agreement with the results read off from Fig. 13 of Ref. [8] 
is quite satisfactorily, even though the statistical errors are obviously already 
quite large. For and some entries in Table 4 are no longer reliable 
and only the orders of magnitude should be trusted. Here we have certainly 
reached the limit of the present simulations, in particular for g = 10, and 
just as a very rough estimate we hnally quote ~ 10^^, 10^®, and 10^®, 
and ~ 10^^, 10^^, and 10^®, for q = 10, 15, and 20, respectively. As an 
example we show in Fig. 2 for g = 10 the cumulant expansion around f3t of the 
energy density e = —{d/dP){—Pf) in the disordered phase and compare the 
results with extrapolations obtained by the standard reweighting method. 

The Table 5 collects the expectation values for the susceptibility Xd at 
/5t, using the two different cluster estimators (12). For comparison we have 
also included the integral over the zero-momentum correlation function g{i) 
[16, 17], which was also computed by employing a Swendsen-Wang cluster es¬ 
timator. One can show that this amounts only to a different implementation 
of precisely the same operations needed to compute directly the Swendsen- 
Wang cluster estimator in eq. (12), and therefore the hrst and second lines 
for each lattice size in Table 5 in fact turn out to be identical as they should. 

The results in Table 6 for the magnetization in the ordered phase clearly 
conhrm that m = m' with high precision. About 5-6 signihcant digits of the 
numerical estimates agree within their la error bounds with the exact results 
of Ref. [12], which provides further evidence that all measured numbers can 
be interpreted as pure phase expectation values. Furthermore we note that 
as expected the higher moments of ilAnax and M^us do not agree. More 
quantitatively we find consistently that Xo > Xo ^'^d The 

proper susceptibility Xo for g = 10 has also been considered in the low- 
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temperature series analysis of Briggs et al. (BEG) [20]. They obtained an 
estimate of x(/3+)beg = 2.44(9), leading to Xo = (^)^x(A’^)beg = 3.01(12), 
which is signihcantly smaller than our value of 4.744(42) on the 300 x 300 
lattice. 


4 Summary 

We have performed high-precision Monte Carlo simulations in the ordered 
and disordered phase of 2D g-state Potts models with q = 10, 15, and 20 
at their hrst-order transition point, working with large lattices of linear size 
L ^ (= 300, 120, 80). As an important self-consistency test the hrst 

three energy cumulants are fonnd to satisfy the duality relations with high 
precision, and both the energy and the magnetization are fully consistent with 
Baxter’s exact values. As our main result we obtain for the hrst three energy 
cnmnlants very good agreement with recent resnmmed large-g expansions 
of Bhattacharya et al. [8], indicating that their techniqne can give reliable 
results at least down to g = 10. Also the fourth- to sixth-order cumulants are 
fonnd in reasonably good agreement, albeit the accnracy of both methods 
decreases with increasing order. We hnd, however, signihcant differences to 
low- and high-temperature series analyses of the specihc heat and magnetic 
susceptibility of the g = 10 Potts model at f3f 
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Appendix 

This appendix lists the relation between cumulants and central moments. 
To simplify the notation we employ here the dehnitions = Vk,^'^ and 
Hn = V= {{E—{E)pY)p. With the help of computer algebra we obtained 

^4 = ^4 3/i2, 

^5 = A'-S ~ 10p2P'3; 

^6 = h-e ~ 15/i2Ai4 — lO/ig + 30/i2, 

^7 = h7 ~ 21/42/^5 ~ 35/43/i4 + 210/i3/i2, 

/tg = Ps ~ 28/42^-6 ~ 56/i3/i5 + 420p4/i2 — 35/44 + 560/43/42 — 630/42, 

/tg = jJ,Q — 36/42/47 — 84/43/4g — 126/44/45 T 756/45/42 

+ 25 20 / 42 / 43/44 + 560/43 ~ 75 60 / 43 / 42 , 

/4l0 = hlO — 45/42/48 — I2O/43/47 — 210/44/4g — 126/4g 

T 1260/4g/42 T 504O/42/43/45 T 315O/44/42 T 42OO/44/43 
— 189 00/44/42 — 378OO/43/42 + 22 6 80/42- 
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Table 1: Integrated autocorrelation time Tint,e of the energy and the number 
of Monte Carlo update sweeps (MCS) in units of rint,e- 



g = 10 
300 X 300 

g = 15 
120 X 120 

g = 20 

80 X 80 

orderec 

. phase (hea 

t-bath algorithm) 

^int,e 

MCS/ri„t,e 


^ 20 

640 000 

^ 9 

1280 000 

disorderec 

. phase (single-cluster algorithm) 

MCS/ri„t,e 

^ 59 

600000 

^ 18 

9 000 000 

^ 25 

4 200 000 


Table 2: Comparison of numerical and analytical results for energy cumulants 
at Pt in the disordered phase. 


Observable 

g = 10 

L = 150 

g = 15 

L = 60 

g = 20 

L = 40 

Cd {d, Lx L) 

-0.96812(15) 

-0.75053(13) 

-0.62648(20) 

Cd {d, 2L X L) 

-0.968190(81) 

-0.750510(65) 

-0.626555(97) 

ed (d, 2L X 2L) 

-0.968186(18) 

-0.7504949(73) 

-0.626519(13) 

Cd (exact) 

-0.968203... 

-0.750492... 

-0.626529... 

Cd (d, L X L) 

18.33(17) 

8.695(47) 

6.144(43) 

Cd {d, 2L X L) 

18.34(12) 

8.665(29) 

6.140(27) 

Cd (d, 2L X 2L) 

18.437(40) 

8.6507(57) 

6.1327(38) 

Cd (large g) 

18.43(2) 

8.657(3) 

6.1326(4) 

Lx L) 

-2010(100) 

-171.0(5.1) 

-54.7(1.9) 

(d, 2L X L) 

-2031(73) 

-176.1(3.8) 

-53.9(1.5) 

(d, 2L X 2L) 

-2015(26) 

-176.01(76) 

-54.85(29) 

(large g) 

-1834(200) 

-174(4) 

-54.7(4) 
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Table 3: Comparison of numerical and analytical results for energy cumulants 
at I3t in the ordered phase. 


Observable 

g = 10 

g = 15 

g = 20 


L = 150 

L = 60 

L = 40 

Co (o, Lx L) 

-1.664177(81) 

-1.765850(34) 

-1.820722(43) 

Co (o, 2L X L) 

-1.664262(57) 

-1.765875(26) 

-1.820689(14) 

Co (o, 2L X 2L) 

-1.664224(58) 

-1.765914(27) 

-1.820659(20) 

Co (exact) 

-1.664253... 

-1.765906... 

-1.820684... 

Co (o, Lx L) 

17.95(13) 

8.016(21) 

5.351(15) 

Co {d, Lx L) 

17.88(17) 

8.037(47) 

5.373(43) 

Co (o, 2L X L) 

17.81(10) 

8.004(19) 

5.3612(55) 

Co (d, 2L X L) 

17.89(12) 

8.007(29) 

5.369(27) 

Co (o, 2L X 2L) 

18.00(10) 

7.990(19) 

5.3608(88) 

Co (d, 2L X 2L) 

17.989(40) 

7.9931(57) 

5.3613(38) 

Co (large q) 

17.98(2) 

7.999(3) 

5.3612(4) 

(o, Lx L) 

1979(87) 

180.5(3.1) 

57.0(1.3) 

/i(3)(d, Lx L) 

2026(100) 

175.7(5.1) 

56.9(1.9) 

fi^^Ho,2Lx L) 

1836(71) 

189.7(5.1) 

56.24(40) 

/i(3) (d, 2L X L) 

2047(73) 

180.8(3.8) 

56.1(1.5) 

/i(3) (o, 2L X 2L) 

2030(110) 

177.2(3.2) 

56.83(71) 

/i(3) (d, 2L X 2L) 

2031(26) 

180.67(76) 

57.09(29) 

(large q) 

1900(200) 

179(4) 

56.9(4) 
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Table 4: Numerical estimates of higher-order energy cumulants in the disor¬ 
dered and ordered phase at (5t. 


Observable 

q = 10 

L = 150 

g = 15 

L = 60 

g = 20 

L = 40 

K'*’{d,2Lx L) 
4"*’ (d, 2L X 2L) 
kW (o, 2L X L) 
kW (o, 2L X 2L) 
(large q) 

1.61(14) X 10® 
1.583(64) X 10® 
1.10(13) X 10® 
1.55(22) X 10® 
1.3(2) X 10® 

3.07(24) X 10^ 
3.193(48) X 10^ 
3.95(47) X 10^ 
2.93(20) X 10^ 
3.1(2) X 10^ 

4.87(38) X 10® 
4.905(89) X 10® 
4.67(14) X 10® 
4.79(22) X 10® 
5.0(1) X 10® 

{d, 2L X L) 
(d, 2L X 2L) 
4®) (o, 2L X L) 
4^) (o, 2L X 2L) 
(large q) 

-2.39(39) X 10^ 
-2.40(26) X 10^ 
1.03(25) X 10^ 
1.98(51) X 10^ 
10® 

-1.10(23) X 10^ 
-1.170(59) X 10^ 
1.73(49) X 10^ 
0.92(15) X 10^ 
^ lO’^ 

-8.1(1.!) X 10® 
-8.42(32) X 10® 
7.16(65) X 10® 
7.37(70) X 10® 

^ 10® 

Kf (d, 2L X L) 
nf (d, 2L X 2L) 
4®) (o, 2L X L) 
(o, 2L X 2L) 
(large q) 

5.1(1.2) X 10^2 
5.8(1.3) X 10^2 
1.13(48) X 10^2 
2.8(1.4) X 10^2 
^ 10^2 

6.6(2.7) X 10® 
7.2(1.1) X 10® 
11.9(5.4) X 10® 
4.5(1.3) X 10® 

3 X 10® 

1.77(34) X 10® 
2.33(16) X 10® 
1.63(36) X 10® 
1.38(26) X 10® 

^ 2 X 10® 

4"^(d, 2Lx L) 
K^P (d, 2L X 2L) 
kP (o, 2L X L) 
kP (o, 2L X 2L) 

-1.20(39) X 10^® 
-2.01(66) X 10^® 
1.15(92) X 10^® 
3.9(4.2) X 10^® 

-5.5(3.5) X 10^2 
-6.8(2.1) X 10^2 
1.03(60) X 10^® 
3.0(1.2) X 10^2 

-4.0(1.!) X 10®° 
-8.47(96) X 10®® 
5.0(2.2) X 10®° 
2.6(1.0) X 10®° 

kP (d, 2L X L) 
kP (d, 2L X 2L) 
kP (o, 2L X L) 
kP (o, 2L X 2L) 

2.9(1.3) X 10^® 
8.6(3.5) X 10^® 

1 X 10^’^ 

^ -2 X 10^® 

5.5(4.4) X 10^® 
8.8(4.4) X 10^® 
1.01(69) X 10^® 
2.1(1.!) X 10^® 

7.9(3.7) X 10®2 
3.57(56) X 10®® 
2.0(1.3) X 10®® 

3.1(4.2) X 10®2 


14 





























Table 5: The magnetic susceptibility Xd at jSt in the disordered phase using 
different estimators. 


Lattice 

Observable 

g = 10 

L = 150 

g = 15 

L = 60 

g = 20 

L = 40 

Lx L 

(g_\)2 Et=i9{^) 
^(|C'n<i,sw/(|C'|)d,sw 
C* )d,sc 

4.224(16) 

4.224(16) 

4.224(16) 

0.7306(14) 

0.7306(14) 

0.7310(14) 

0.3092(58) 

0.3092(58) 

0.3090(58) 

2L X L 

C'Hd.Sw/ ( C* )d,SW 

^(|C'|)d,sc 

4.2306(89) 

4.2306(89) 

4.2327(89) 

0.73093(68) 

0.73093(68) 

0.73094(65) 

0.30954(32) 

0.30954(32) 

0.30952(31) 

2L X 2L 

;7^(|CP)d,sw/( c ).i,sw 

4.2326(18) 

4.2326(18) 

0.730386(79) 

0.730386(79) 

0.309356(39) 

0.309356(39) 
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Table 6: The magnetization m and m' at (3t in the ordered phase, using the 
two estimators M^ax and Mdus, and the corresponding susceptibilities and 
third moments. 


Observable 

g = 10 

g = 15 

g = 20 


L = 150 

L = 60 

L = 40 

m { L X L) 

0.857047(71) 

0.916631(21) 

0.941199(21) 

m' {Lx L) 

0.857047(71) 

0.916634(21) 

0.941197(21) 

m {2L X L) 

0.857113(49) 

0.916648(16) 

0.9411782(66) 

m' {2L X L) 

0.857113(49) 

0.916648(16) 

0.9411791(66) 

m {2L X 2L) 

0.857081(49) 

0.916672(15) 

0.9411694(97) 

m' {2L X 2L) 

0.857077(49) 

0.916672(15) 

0.9411694(97) 

m (exact) 

0.857106... 

0.916663... 

0.9411759... 

Xo (Lx L) 

4.750(60) 

0.8090(36) 

0.3348(17) 

X'„ {Lx L) 

4.704(60) 

0.7989(36) 

0.3305(17) 

X 

o 

do 

X 

4.663(43) 

0.8095(38) 

0.33509(55) 

(2ix L) 

4.623(43) 

0.7997(38) 

0.33076(55) 

Xo {2L X 2L') 

4.744(42) 

0.8052(28) 

0.33551(81) 

X'o {2 Lx2L) 

4.700(43) 

0.7953(28) 

0.33118(80) 

mf'> {Lx L) 

-1521(85) 

-45.9(1.2) 

-8.55(32) 

mf\ Lx L) 

-1505(84) 

-45.4(1.2) 

-8.44(32) 

{2L X L) 

-1372(62) 

-49.4(2.2) 

-8.321(88) 

mf\2Lx L) 

-1362(62) 

-48.9(2.2) 

-8.216(88) 

{2L X 2L) 

-1532(74) 

-45.3(1.0) 

-8.44(13) 

mf\2L X 2L) 

-1517(75) 

-44.8(1.0) 

-8.33(13) 
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Figure 1: The energy probability distribution P{e) of the 10-state model at 
Pt in the ordered and disordered phase for L x L (dotted lines), 2Lx L (solid 
lines), and 2L x 2L (dashed lines) lattices. The area under each peak is 
normalized to unity. 



P 

Figure 2: Cumulant expansion of the energy density e = —{d/dP){—Pf) 
for g = 10 around Pt = log(l -|- vTO) ^ 1.426 062... in the disordered 
phase. The numbers n = 3,4, 5,... in the legend indicate the highest order 
of the cumulants involved in the expansion. The open circles show for 
comparison the energy as obtained by standard reweighting. 
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